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The laser-driven acceleration of high quality proton beams from a double-layer target, comprised 
of a high-Z ion layer and a thin disk of hydrogen, is investigated with three-dimensional particle- 
in-cell simulations in the case of oblique incidence of a laser pulse. It is shown that the proton 
beam energy reaches its maximum at a certain incidence angle of the laser pulse, where it can be 
much greater than the energy at normal incidence. The proton beam propagates at some angle with 
respect to the target surface normal, as determined by the proton energy and the incidence angle. 
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The method of charged particle acceleration by using 
laser light is very attractive, since the acceleration rate 
is much higher and the facility size can be substantially 
smaller than standard accelerators. Ion acceleration dur- 
ing the high intensity electromagnetic wave interaction 
with plasmas was proposed more than 50 years ago pQ. 
Currently, ion acceleration experiments using high power 
lasers close to petawatt levels are going on all over the 
world [2 j. Laser driven fast ions are considered in regard 
to applications ranging from hadron therapy [3], fast ig- 
nition of thermonuclear targets [4], production of PET 
sources [5], conversion of radioactive waste [6 j, a laser- 
driven heavy ion collider [7], injectors for standard accel- 
erators [8], and proton radiography [9 to proton dump 
facilities for neutrino oscillation studies [10] (see Refs. 
[TT] and literature quoted therein). 

The typical energy spectrum of laser accelerated par- 
ticles from unoptimized targets is thermal-like, with a 
cut-off at a maximum energy. On the other hand, al- 
most all the above mentioned applications require high 
quality proton beams, i.e. beams with sufficiently small 
energy spread A£/£ < 1. As suggested in Ref. [12], 
such a required beam of laser accelerated ions can be 
obtained using a double-layer target, which consists of 
high-Z atoms and a thin coating of low-Z atoms. Exten- 
sive computer simulations of this target were performed 
in Refs. 13 and [14 , where multi-parametric particle- 
in-cell (PIC) simulations were used to optimize the laser- 
driven proton acceleration by choosing appropriate laser 
and target conditions. The feasibility of the double-layer 
target scheme was demonstrated experimentally with mi- 
crostructured targets in Ref. [15 . The effects of target 
shaping on the laser-driven ion acceleration were also re- 
ported in Refs. [16 . Previously, the double-layer tar- 
get scheme for high and controllable quality ion accel- 
eration has been mostly studied in the configuration of 
normal incidence of the laser pulse on the target. How- 
ever, the case of oblique incidence provides an additional 
parameter for manipulation of the fast ion energy, the 
emittance, energy spectrum and the proton beam prop- 
agation direction. As is well known, the energy transfer 



from a p-polarized obliquely incidence electromagnetic 
wave to the electron energy via the so-called "vacuum 
heating" mechanism [17 depends on the incidence an- 
gle and is substantially higher than for normal incidence. 
A stronger electron heating results in a stronger electric 
field generation due to the electric charge separation ef- 
fect, which in turn leads to more efficient ion acceleration. 

In this Letter, we study the dependence of the ion 
beam energy and quality on the laser incidence angle. We 
use an idealized model, in which a Gaussian p-polarized 
laser pulse is incident on a double-layer target of colli- 
sionless plasmas. The simulations are performed with 
a three-dimensional massively parallel electro magnetic 
code, based on the PIC method [18]. The number of grid 
cells is equal to 2800 x 720 x 720 along the x, y, and z 
axes. The total number of quasi-particles equals 7 x 10 7 . 
The size of the simulation box is 100A x 25. 5A x 25. 5A, 
where A is the laser wavelength. The boundary condi- 
tions for the particles and for the fields are periodic in the 
transverse (y,z) direction and absorbing at the bound- 
aries of the computation box along the x axis. Here the 
laser wavelength determines the transformation from di- 
mensionless to dimensional quantities and vice versa. Be- 
low, the dimensional quantities are given for A = 0.8/im; 
the spatial coordinates are normalized to A and the time 
is measured in the laser period, 2tt/uj. 

The Gaussian laser pulse with the dimensionless ampli- 
tude a = eEo/m e coc = 30, which corresponds to the laser 
peak intensity 2 x 10 21 W/cra 2 , is 8A long in the propa- 
gation direction and it is focused to a spot with size 6A 
(FWHM). The oblique incidence of the laser pulse is re- 
alized by tilting the target around the z— axis, Fig. [ljd) , 
while the laser pulse propagates along the x— axis. Both 
layers of the double-layer target are shaped as disks. The 
first, gold, layer has a diameter 10A and thickness 0.5A. 
The second, hydrogen, layer is narrower and thinner; its 
diameter is 5A and thickness is 0.03A. The electron den- 
sity inside the gold layer is n e = 1.6 x 10 22 cm -3 ; inside 
the proton layer it is n e = 5 x 10 20 cm -3 . 

Figure [I] shows the proton beam acceleration for two 
cases of normal and oblique incidence. For the present 
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FIG. 1: Proton acceleration driven by the laser pulse with 
normal (a,b) and oblique (c,d) incidence. (a,c): Electric field 
magnitude (isosurface for value a = 2) at t=10; half of the 
box is removed to reveal the internal structure. (b,d): Dis- 
tribution of gold ions (blue), electrons (yellow), and protons 
(light blue) at t=80. The laser pulse incidence angle in (c,d) 
is6> = 30°. 



simulation parameters, the target is partially transpar- 
ent for the electromagnetic wave (due to the relativistic 
transparency effect), which, according to Refs. [13] and 
[T4] , corresponds to an optimal condition of the ion accel- 
eration. In the case of oblique incidence, a substantially 
larger portion of the electrons is blown off the target due 
to the effect of "vacuum heating." 

In order to examine the dependence of the energy 
achieved by the protons on the incidence angle, we 
performed extensive 2D and 3D PIC simulations. Fig. [2] 
shows the proton energy as a function of the laser pulse 
incidence angle, normalized by the maximum achievable 
energy. We see that the maximum proton energy is 
reached at the incidence angle approximately equal to 
30°, for which case we undertook 3D PIC simulation 
shown in Fig. [TJc,d). The proton energy value is ap- 
proximately doubled at #=30° compared to the case of 
normal incidence, #=0°. This also can be seen in the in- 
set in Fig. [2] showing the proton energy spectra at t=80. 
While the proton energy increased from ~20 MeV at nor- 
mal incidence to ~45 MeV at oblique incidence with an 
optimal angle, the energy spread is also increased from 
10% to 23%. 

If we invoke the above mentioned "vacuum heating" 
mechanism for the laser pulse energy transformation into 
the electron component energy, and then, via acceleration 
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FIG. 2: Proton energy, normalized by the maximum, vs the 
laser pulse incidence angle. In the inset: proton energy spec- 
tra at t=80, as obtained in 3D PIC simulation. 



in the electric field of the charge separation, into the en- 
ergy of fast protons, we find that the laser pulse energy 
deposited to the target and the efficiency of pushing the 
electrons out from the target depend on the incidence 
angle 0. The count erplay between these two effects leads 
to the formation of the fast proton energy maximum at a 
certain incidence angle. Under the conditions of our sim- 
ulations this angle is approximately equal to 30°. The 
analysis of the time evolution of the proton energy spec- 
trum shows that both the average energy and energy 
spread increase with time. The energy spread appears 
to grow in the oblique incidence case, which can be ac- 
counted for by the asymmetry of the quasistatic electric 
field along the target surface at oblique incidence. We 
note that, according to Ref. [12], the energy spread can 
be decreased by reducing the thickness and transverse 
size of the low-Z ion (proton) layer. 

Since the proton acceleration occurs due to the elec- 
tric field generation through the electron ejection from 
the target, it is important to analyze the dependence 
of the number of the electrons leaving the target under 
the laser radiation action on the incidence angle. Fig- 
ure [3] shows the ratio of electrons pushed out from the 
target versus the incidence angle at time t=20. The 
dependence at t=20 is chosen because at this time the 
number of ejected electrons saturates as seen in the in- 
set in Fig. [3] The ratio of ejected electrons as a func- 
tion of the incidence angle 6 is defined by the expression 



Ve = (N e (0)-N eirnin )/(N e 



i), where N e (0) is 



the number of ejected electrons, N eiTnax = max{iV e (#)} 
and Ne^min = min{7Ve(#)} are the maximum and mini- 
mum values, respectively. We see that this ratio reaches 
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FIG. 3: Escaping electron ratio r\ e vs the laser pulse incidence 
angle at t=20. In the inset: normalized number of electrons 
swept off the target vs time, as obtained in 3D PIC simulation. 



its maximum at the incidence angle about 30°, which ev- 
idences a strong correlation between the dependences of 
the number of ejected electrons (Fig. [3| and the pro- 
ton energy (Fig. [2| on the laser pulse incidence an- 
gle. We see from Fig. [3] that the efficiency of the 
electron ejection is substantially higher for #=30° than 
in the case of normal incidence. As a consequence, a 
stronger electric field is produced and the protons are 
accelerated to higher energy. This is also illustrated 
in Fig. [4j by the correlation between the average pro- 
ton energy and the ejected electron number. Here the 
proton energy ratio at 9 is defined by the expression 
Vs 



is the average proton energy, E p ^ r 



p,rmn 



S Pi min), where £ p (0) 
x = max{£p(#)} and 
min{£ p (#)} are the maximum and minimum 
values, respectively. We see that the fast proton energy is 
roughly proportional to the number of ejected electrons. 

As seen from Fig. [T] the accelerated proton bunch re- 
tains the form of the disk. In the case of the normal 
incidence, the proton disk moves in the direction normal 
to the target surface while its surface remains parallel 
to the first (high-Z ion) layer. In the case of oblique in- 
cidence, the proton disk motion direction deflects from 
the target normal by a noticeable angle, 0, while the 
disk surface is tilted with respect to the high-Z ion layer, 
Fig. [ljd) . We note that a deflection of the accelerated 
proton bunch has been observed in the experiments on 
the laser - solid target interaction [19,. Under the con- 
ditions of our simulations, which correspond to the rela- 
tively higher laser intensity, the deflection and tilt can be 
explained by relativist ic effects. As is known, the Lorentz 
transformation to the frame of reference moving with the 
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FIG. 4: Proton energy ratio vs escaping electron ratio at 
t=20. 



velocity V = — csin# along the target surface, i.e. the 
gamma factor is given by 7 = 1/ y/l — V 2 /c 2 = 1/ cos#, 
changes the configuration of the wave-target interaction 
from oblique to normal incidence [20 , so that the wave 
frequency and wave vector of the incident electromag- 
netic wave become oj' = uo cos 6 and k' = fc cos 0e|| , where 
e|| is the unit vector along k 7 . In this new boosted refer- 
ence frame, the protons have a transverse component of 
momentum equal to p' ± = —m p c tan where m p is the 
proton mass. As a result of the acceleration in the elec- 
tric charge separation field, the protons acquire the lon- 
gitudinal momentum pj. and their energy becomes equal 

to £p = yjm 2 c 4 + p f2 c 2 + p f2 c 2 . Performing the Lorentz 

transformation back to the laboratory reference frame, 
we obtain that p\\ = pj,, p±_ = j(p f u — £' p V/c 2 ) and the 
deflection angle (j) = arctan(p^/p||) is equal to 



arctan 



(i) 

In the limit p\\/m p c <C 1, which corresponds to the pa- 
rameters of our simulations, this expression yields <fi « 
^£ p /2m p c 2 sin#. For p\\/m p c ^> 1, the angle (j) becomes 
(j) ~ i.e. the protons are accelerated almost along the 
laser pulse propagation direction. 

In order to account for the proton disk tilting, we note 
that the obliquely incident laser pulse front propagates 
along the target surface with a superluminal velocity 
Vf = cj sin#. The time delay between moments when 
different disk elements separated by the distance A/ start 
to move is equal to At = Al smO/c. The displacement of 
the proton layer elements in the direction of the target 
normal can be estimated as A£|| = p\\At/m p . This gives 
the angle of the proton disc tilt, x — arctan(A£||/AZ), i.e. 
X ~ y/ 28 p / m p c 2 sin 0. This effect is seen in the simula- 
tions at an early time of the proton acceleration. Later, 
higher dimensional effects come into play and the tilting 
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FIG. 5: Deflection angle of the proton bunch vs the laser 
pulse incidence angle at t=80. The solid curve is given from 
Eq. 0. 

angle changes. 

The dependence of the angle of the deflection of the 
proton motion from the target normal, 0, on the laser 
pulse incidence angle, is seen in Fig. [5] Here the 
deflection angle, </>, is defined as the angle between the 
normal to the target surface placed at the target center 
and the average radius-vector from the target center to 
the proton layer, which equals S^x^/A/^, where is 
the radius- vector of the z-th proton out of N p protons in 
the accelerated beam. As seen from Eq. ([!]), in the limit 
of a small incidence angle, the proton bunch is accel- 
erated almost along the direction normal to the target 
with <j> being a linear function of 6. When the incidence 
angle, increases, the growth of the deflection angle, (/>, 
saturates, in agreement with the results of simulations, 
shown in Fig. [2] 

In conclusion, we have found that the proton accel- 
eration during laser pulse interaction with double-layer 
targets is more efficient for oblique incidence than for 
normal incidence. It is shown that the proton beam en- 
ergy reaches its maximum at a certain incidence angle of 
the laser pulse, where it can be much greater than the en- 
ergy at normal incidence. The proton beam propagates 
at some angle with respect to the target surface normal, 
as determined by the proton energy and the incidence 
angle. In the limit of nonrelativistic proton energy, the 
deflection angle is relatively small. However, its value 
« 6° (see Fig. |5| results in the proton beam deflection of 
the order of 10 cm after the protons have propagated over 
1 m distance. Therefore this angle should be taken into 
account in the real facility design to make proper mutual 
location of the proton source and the target. This pro- 
vides a way to control the proton energy and the direction 
of the proton beam propagation by adjusting the inci- 
dence angle of the laser pulse. Such flexibility and speed 
in controlling the proton directionality and energy will be 
invaluable for conformal cancer therapy treatments such 



as pixel scanning using heavy ions [21]. 
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